Accepted for publication in ApJ 

Preprint typeset using I^T^X style cmulatcapj v. 04/20/08 



OS 
O 
O 

> 

O 

O 

in 



< 

6 
to 



O 



o> 
o 



X 



EMBEDDED YOUNG STELLAR OBJECT CANDIDATES IN THE ACTIVE STAR FORMING COMPLEX W51: 

MASS FUNCTION AND SPATIAL DISTRIBUTION 

Muu KANG 1 ' 2 ' 3 , John H. Bieging 3 , Matthew S. Povich 4 ' 5 , and Youngung Lee 1 

Accepted for publication in ApJ 

ABSTRACT 

We present 737 candidate Young Stellar Objects (YSOs) near the W51 Giant Molecular Cloud 
(GMC) over an area of 1.25° x 1.00° selected from Spitzer Space Telescope data. We use spectral 
energy distribution (SED) fits to identify YSOs and distinguish them from main-sequence or red giant 
stars, asymptotic giant branch stars, and background galaxies. Based on extinction of each YSO, we 
separate a total of 437 YSOs associated with the W51 region from the possible foreground sources. We 
identify 69 highly embedded Stage 0/1 candidate YSOs in our field with masses > 5 M Q (corresponding 
to mid- to early-B main-sequence spectral types) 46 of which are located in the central active star 
forming regions of W51A and W51B. From the YSOs associated with W51, we find evidence for mass 
segregation showing that the most massive YSOs are concentrated on the W51 H II region complex. 
We find a variation in the spatial distribution of the mass function (MF) of YSOs in the mass range 
between 5 Mq and 18 Mq. The derived slopes of the MF are —1.26 and —2.36 in the active star- 
forming region and the outer region, respectively. The variation of the MF for YSOs embedded in 
the molecular cloud implies that the distribution of stellar masses in clusters depends on the local 
conditions in the parent molecular cloud. 

Subject headings: H II regions — stars: formation — ISM: individual (W51) — infrared: ISM 



1. INTRODUCTION 

Stars form in various environments. High-mass stars 
are especially important because they affect their en- 
vironment through such phenomena as outflows, stellar 
winds, and supernovae. To understand the interaction of 
stars and the interstellar medium (ISM) it is necessary 
to identify embedded young stellar objects (YSOs) in the 
parent molecular cloud. Therefore taking a complete cen- 
sus of YSOs is an important step toward understanding 
YSOs and their environment on a large scale. 

W51 is an active star-forming region located in the 
Sagittarius spiral arm. The distance is uncertain, 
with published values including 5.lj^4 kpc (|Xu et al.1 
12001 . 6.1 ± 1.3 kpc (llmai et al.1 120021 or 8.3 ± 2.5 
kpc ( Schneps et al.lll981f ). The radio continuum sources 
comprising W51 are within a massive giant molecu- 
lar cloud (GMC) (jMufson fc Liszti Il979f l. Observations 
showing the star forming activity have been carried 
out in various wavelengths of radio continuum, HI line, 
mole c ular lines, water masers, n e ar-IR, X-rav dBiegineJ 
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sive star formation 
and therefore a good place to study the feedback between 
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stars and the ISM. 

To understand the interaction of molecular clouds and 
YSOs in the W51 H II region complex we mapped a 1.25° 
x 1.0° region in the J = 2 — 1 transitions of 12 CO and 
13 CO with the 10 m Heinrich Hertz Telescope (HHT) 
on Mount Graham, Arizona. The molecular line maps 
will be presented in a separate paper (Bieging et al., in 
preparation). In the present paper, we use Spitzer data 
for identifying the YSO candidates associated with the 
molecular clouds. Spitzer surveys in the mid-IR allow 
sensitive studies of deeply embedded star formation. Al- 
though extinction is significantly reduced in the mid-IR, 
it is still important for active star forming regions in 
giant molecular clouds. In particular, the W51 region 
shows no detecta ble increase in star co unts because of 
severe extinction (| Benjamin eTal1l2TJ05lL Therefore wo 
have to consider extinction effects carefully in identify- 
ing and characterizing candidate YSOs. Combining the 
YSOs selected from the Spitzer photometry and molec- 
ular clouds with kinematic information, we can examine 
the feedback process between star(s) and ISM. The de- 
tailed analysis for the interaction of YSOs and associated 
molecular clouds w ill be presented in other publications 
(|Kang et al.l (|2009D : Kang et al., in preparation). In this 
paper, we tabulate all of the YSO candidates and discuss 
their spatial distribution and mass function. In Section[2l 
we introduce the data sets used and describe our meth- 
ods for selecting and classifying candidate YSOs. We 
discuss some of the implications of our results in Section 
[3] and summarize the results in Section [4] 

2. YSO SELECTION AND CLASSIFICATION 
2.1. Data 

The Galactic Legacy In frared Mid-Plane Sur vey Ex- 
traordinaire (GLIMPSE I; iBeniamin et al.ll2003f) survey 
covered the Galactic plane (10° < \l\ < 65°, |6| < 1°) 
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TABLE 1 
Source Counts in the W51 Region 3 



Sources 



Number 



Fig. 1. — A color image of the W51 complex composed of MIPS 
24 /im (red) and IRAC 8.0 (im (green). Hot dust grains and PAH 
dominated emission appear in red and green, respectively. The 
contours show 21 cm radio continuum emission ma de by combin- 
ing V LA data and Effelsberg 100-m telescope data I IKoo Ik Moonl 
HWD . The contour levels are -0.015, 0.015, 0.05, 0.10, 0.3, 0.5, 
1.0, 1.5, and 2.4 Jy bea m -1 . Cross symbols mark the radio con- 
tinuum sources listed by[Koo & Moon ( 19 971) . 



with the four mid-IR bands (3.6, 4. 5, 5.8, and 8.0 nm ) of 
the Infrared Array Camera fIRAC: lFazio et al.ll2004[ ) on 
the Spitzer Space Telescope. Each IRAC band contains 
different spectral features : Band 1 [3.6] shows mainly 
continuum emission from stars, Band 2 [4.5] traces H2 
rotational transitions arising in the shocked gas associ- 
ated with outflows, and emission in Bands 3 [5.8] and 
4 [8.0] is dominated by polycyclic aromatic hydrocarbon 
(PAH) features. For this study, we have retrieved images 
of the 1.25° x 1.0° region of W51 centered on (l,b) = 
(49.375°, —0.2°) (note that this is the region we mapped 
in CO) by combining two mosaic images which have a 
high resolution (l'/2 pixel" 1 ) in all four IRAC bands (Fig- 
ure [TJ. For this same region, we used the GLIMPSE I 
Catalog including the sources detected at least twice in 
one band with a S/N > 5 and at least once with a S/N 
> 3 in an adjacent band. Generally, the flux density lim- 
its insure that the detections are > 5a. The 5er point- 
source detection limits of the GLIMPSE are nominally 
0.2, 0.2, 0.4, and 0.4 mJy for the IRAC 3.6, 4.5, 5.8, 
and 8.0 fim bands, respectively. These limits are signif- 
icantly higher in regions of bright diffuse emission. The 
GLIMPSE I Catalog also tabulates JH K R flux densities 
from the 2MASS point source catalog (Skrutski e et al.l 
[20061) for all GLIMPSE sources with 2MASS identifica- 
tions. 

MIPSGAL (jCarev et al.ll2005fl is a legacy program cov- 
ering the inner Galactic plane, 10° < \l\ < 65° for |£>| < 
1°, at 24 and 70 /jm with the Multiband I maging Pho- 
tome ter for Spitzer Space Telescope (MIPS; iRieke et al~l 
2004). The resolution of the 24" jjxa mosaics from the 
MIPSGAL survey is 2'.' 4 pixel" x . Wc extracted 24 ^m 
point-sources with F/5F > 7, then bandmerged the 24 
/zm sources with the GLIMPSE Catalog sources using a 
2'.'0 correlation radius. The final source list consists of all 
8 bands combined: 2MASS, GLIMPSE, and MIPSGAL 
24 /xm. A total of 104,582 sources within the 1.25° x 1.0° 
region of W51 centered on (/,&) = (49.375°, -0.2°) were 
selected from the GLIMPSE I Catalog. 



In GLIMPSE Catalog 104,582 

Fit with SED models (N data > 4) 70,767 

Well-fit by stellar photosphere SEDs (x 2 /N data < 4) 69,316 

Poorly-fit by stellar photosphere SEDs (x 2 /N data > 4) 1,451 

Well-fit by YSO SEDs { X 2 /N data < 4) 1,327 

In final sample of YSO candidates 737 

a The 1.25°x 1.0°rcgion centered on (l,b) = (49.375°, -0.2°) 

2.2. Color Selection 

The colors of most main sequence (MS) stars are near 
zero for any combination of the IRAC and MIPS bands, 
while YSOs have red colors due to the emission from 
surrounding warm dust. As an initial appr oach, we ap- 
lied a well defined color selection criteria of lSimon et al.l 
2007f ) to our source list for identifying YSO candidates 
associated with the molecular cloud around W51. After 
selecting candid ate YSOs from the region in color-color 
space defined bv lSimon et al.l (|2007f ). we classified them 
with the spectral index a in the IRAC bands. The spec- 
tral index is 

d log [AF(A)] 
a logA 

where A is the wavelength and F{\) is the flux at that 
wavelength. We obtain a from a linear fit to the loga- 
rithm of the fluxes for all 4 IRAC bands and then clas- 
sify YSOs as Class I sources with 0.3 < amAc, " Flat 
" sources with —0.3 < Chirac < 0.3, Class II sources 
with — 1.6 < olt rac < —0-3, and Class III sou rces with 
aiRAC < -1.6 (|Ladalll987tlGreene et aD 19941 ). 

Figure ^ a) shows the distribution of YSO s within the 
color selection criteria of lSimon et alJ (|2007f ) represented 
by the solid line on the IRAC [3.6] - [4.5] vs. [4.5] - [8.0] 
color-color diagram. In Figure^ a), a large number of 
YSO can didates with red c olors were not selected be- 
cause the lSimon et al.l (|2007f ) color criteria were designed 
to identify a relatively clean sample of YSOs. Many red- 
dened main sequence stars in our W 51 field were identi- 
fied as YSOs, however, because the iSimon et al.1 (|2007l ) 
color criteria were derived from the Small Magellanic 
Cloud which is not affected by the large interstellar ex- 
tinction of the Galactic plane. The W51 region is sev- 
eral kpc from the Sun in the Galactic plane so inter- 
stellar extinction is significant. Figure HJ&) shows the 
IRAC [3.6] -J4.5] v s. [5.8] - [8.0] color-color diagram by 
ISimon et al.1 ([2007) with YSO categories based on a for 
comparing with the distribution of YSOs found by other 
selection criteria in the same color-color space. 

Figure c ) shows the distri bution of YSOs identifie d 
with the IRAC color criteria of IGutermuth et "all (2008) . 
They classified the YSO candidates into two groups af- 
ter removing sources dominated by PAH emission, shock 
emission and broad- line AGNs. One group is Class I pro- 
tostars and the other is the mor e evolved Class II sources . 
Applying color classification of IGutermuth et aLl ([2008), 
we found that the YSOs candidates in the W51 re- 
gion separate into 366 Class I sources and 1457 Class II 
sources. Although these color criteria selected the inter- 
esting YSO candidates having very red colors, the sample 
also contained a large number of reddened stellar photo- 
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Fig. 2. — (a) IRAC [3.6] — [4.5] vs. [4.5] — [8.0] color-color diagram by the color selection criteria of Simon et al. (2007). All sources were 
detected in all 4 IRAC bands in the GLIMPSE Catalog. YSO candidates classified by the spectral index a arc marked in red for Class I 
objects, orange for the flat-s pectrum sour c es, ye llow for the Class II sources, and cyan for the Class III sources. The solid line represents 
th e color select i on cri teria of Simon ct al. (2007). (6) IRAC [3.6] — [4.5] vs. [5.8] — [8.0] color-color diagram by t he color selection criteri a 
of lSimon et~aLl |(200W) . (c) IRAC [3.6] - [4.5] vs. [5.8] - [8.0] color-color diagram by the color selection criteria of lGutermuth et al.l <|2008h . 
Class I objects are marked in red and Class II objects are marked in yellow, (d) IRAC [3.6] — [4.5] vs. [5.8] — [8.0] color-color diagram with 
YSO candidates identified using the SED fitting tool. Small black dots are sources well-fit by stellar photosphere SEDs. YSO candidates 
are marked in red for Stage 0/1, yellow for Stage II, cyan for Stage III, and green for ambiguous sources. AGB candidates are marked by 
blue asterisk symbols. Blue dot s are well-fit YSO candid ates that were discarded from the final sample. Reddening vector for Av=20 mag 
based on the extinction laws of Indcbctouw ct al. (2005) is shown as a filled arrow. Black crosses show typical photometric errors. 



spheres because the color criteria were derived from the 
study of the nearby star forming region, NGC 1333. 

Class I, Flat, and Class II sources in Figure[U&) belong 
to Class I and relatively red Class II samples in Figure 
[2]Jc). Color selection criteria may be suitable as a method 
to select the objects in a relatively young stage. However, 
more evolved sources, which include the Class III objects 
in Figure WLJ>) and those Class II objects around zero 
color in Figure Wi c )i may be mis-identified because of 
the large interstellar extinction toward the W51 region. 



2.3. SED fitting 



We us ed the spectral energy distribution (SED) fitting 
tools of iRobitaille et afl (|2007l ) for identifying and clas- 
sifying YSOs in the W51 giant molecular cloud complex. 
The SED fitting tool works as a regression method to 
find the SEDs within a specified x 2 from a large grid 
of models after fitting the input data points. The grid 
of models contains stellar atmospheres, a limited set of 
spectra of galaxies, AGB stars, and YS O models. The 
grid o f YSO models was computed by IRobitaille et al.1 
(2006) using the 20 , 000 2-D radiat ion transfer models 
from IWhitnev et al.l (|2003allbl . l200l . Each YSO model 
has SEDs for 10 viewing angles (inclinations), so the 
total YSO grid consists of 200,000 SEDs. The SED 
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TABLE 2 
Magnitudes of YSO Candidates 



No. IRAC designation Magnitudes Assoc. a 











u 

ll 




v\ 6 


Id Kl 


K 8 

[o.SJ 


8 n 


\')A\ 




1 


SSTGLMC 


G048.7567-00.6341 


15.89 


14.81 


14.00 


12.58 


12.15 


11.85 


10.94 


7.43 


F 


2 


SSTGLMC 


G048.7579-00.2797 


14.17 


13.37 


12.94 


12.40 


12.30 


11.66 




6.03 


F 


3 


SSTGLMC 


G048.7605-00.0388 








13.05 


12.28 


11.07 


10.02 


6.24 


VV 














14.33 


13.20 


11.68 


10.94 


7.62 


vv 

V V 


5 


SSTGLMC 


G048.7637+00.2022 








13.80 


12.37 


11.27 


10.67 


6.47 


VV 


6 


SSTGLMC 


G048.7655+00.1017 


14.91 


14.08 


13.41 


12.74 


12.64 






8.44 


F 


7 


SSTGLMC 


G048.7667+00.0766 


15.97 


14.10 


12.82 


11.33 


10.80 


10.31 


9.84 


7.60 


F 


8 


SSTGLMC 


G048.7702-00.1505 








13.58 


12.15 


11.46 


11.06 


6.26 


VV 


9 


SSTGLMC 


G048.7703+00.0786 








12.41 


10.88 


9.67 


8.85 


5.81 


w 


10 


SSTGLMC 


G048.7720-00.4792 


16.27 


14.38 


12.90 


10.94 


10.53 


10.03 


9.25 




F 


11 


SSTGLMC 


G048.7726-00.2215 






13.20 


12.06 


11.74 


9.94 




2.28 


F 


12 


SSTGLMC 


G048.7744+00.1091 






14.87 


11.70 


10.47 


9.72 


9.13 


7.20 


W 


13 


SSTGLMC 


G048.7750-00.1507 








13.22 


12.08 


10.41 




5.34 


W 


14 


SSTGLMC 


G048.7772+00.2513 




14.94 


13.88 


12.92 


13.63 


11.11 


9.80 


3.72 


VV 


15 


SSTGLMC 


G048.7780+00.2275 




15.19 


14.29 


13.20 


12.72 


12.15 


11.56 




F 


Note. 


— Table 2 is 


published in its entirety 


in the electronic edition of the Astrophysical Journal. A portion 


is shown 



here for guidance regarding its form and content. 

a Association by the foreground extinction As/ - W, W51 source; F, foreground source 



TABLE 3 

Model parameters for YSO candidates 



No. 


Name 


A v (mag) Ah (M ) 


Ltot (£©) Evol. 


Class c Class d 




(Gl + b) 


L a Ave b U a L a Ave b U a 


L a Ave b U a Stage 


Q I RAC 



1 


G048.7567- 


-00.6341 


2.9 


4.0 


5.4 


2.6 


3.0 


3.1 


13 


74 


86 


11 


11 


II 


2 


G048.7579- 


-00.2797 


0.3 


1.6 


2.6 


1.4 


2.7 


4.0 


17 


33 


45 


0/1 




F 


3 


G048.7605- 


-00.0388 


23.1 


33.3 


42.8 


3.3 


4.1 


4.8 


87 


269 


403 


11 


I 


1 


1 


G048. 76184 


-00.0627 


43.5 


51.8 


60.0 


3.3 


4.0 


4.7 


122 


274 


383 


11 


1 


F 


5 


G048. 76374 


-00.2022 


6.6 


25.7 


53.8 


2.2 


3.5 


4.7 


49 


155 


311 


0/1 


1 


I 


6 


G048. 76554 


-00.1017 


1.6 


3.0 


4.3 


2.8 


3.2 


3.8 


16 


66 


122 


Amb 




11 


7 


G048. 76674 


-00.0766 


7.1 


7.5 


8.1 


3.3 


1.0 


1.5 


129 


229 


327 


II 


11 


11 


8 


G048.7702- 


-00.1505 


1.2 


18.8 


39.6 


1.8 


3.9 


5.7 


35 


106 


188 


0/1 


F 


1 


9 


G048. 77034 


-00.0786 


46.5 


52.7 


60.0 


5.1 


6.6 


7.8 


658 


1788 


2668 


II 


I 


F 


10 


G048.7720- 


-00.4792 


6.2 


9.3 


11.4 


3.8 


1.6 


5.5 


175 


349 


456 


11 


11 


F 


11 


G048.7726- 


-00.2215 


1.1 


7.9 


14.7 


4.8 


6.0 


7.2 


190 


472 


753 


0/1 




I 


12 


G048. 77444 


-00.1091 


24.0 


34.7* 


40.1 


5.4 


5.8* 


6.0 


1179 


1464 


1536 


11 


F 


F 


13 


G048.7750- 


-00.1507 


4.9 


32.2 


56.8 


1.3 


4.3 


6.0 


46 


663 


809 


II 




1 


11 


G048. 77724 


-00.2513 


13.7 


36.8 


58.9 


2.0 


3.6 


5.2 


98 


238 


389 


0/1 


I 


I 


15 


G048. 77804 


-00.2275 


3.7 


8.2 


11.8 


1.4 


2.6 


3.7 


15 


35 


43 


Amb 


11 


11 



Note. — Table 3 is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for 
guidance regarding its form and content. 

a The range of values reported here on the extinctions, masses, and luminosities for each parameter are the lower and upper limit 
of 68% probability for all of the acceptable YSO models. 

Asterisks indicate ranges of 95% probability. 
c YSO class based on spectral slope a using the flux densities detected in all four IRAC bands. 
^ YSO class based on spectral slope a using the flux densities available between 2 and 24 ^m. 



fitting tool provides the evolutionary stage and physi- 
cal parameters such as disk mass, disk accretion rate, 
and stellar temperature of YSOs. The SED fitting tool 
has been used to study st ar forming regions in various 
environments, e.g., M17 dPoyich et al.l 12009). the Ea- 



gle Nebula firTd cbctouw" et al.ll2007D .rRAS 18507+0121 
Shepherd et al.ll2007ft . and the Small Magellanic Cloud 



Simon et al.ll2007h~ 



Table [Tj summarizes the number of sources. We reset 
the photometric uncertainties to a value of 8F/F = 10 % 
for all GLIMPSE Catalog sources, and to 15 % for the 
MIPS 24 fxm fluxes including systematic errors based on 
the GLIMPSE documents 6 . These values were adopted 
instead of the formal errors in the Catalog in order to 
fit without possible bias caused by underestimating the 

6 |http://irsa.ipac.caltech.edu/data/SPITZER/GLIMPSE/doc/| 



flux uncertainties. A total of 70,767 sources for SED fit- 
ting was detected in Ndata > 4 of the 8 IR bands in the 
GLIMPSE Catalog combined with MIPS 24 /im fluxes. 
Our final goal is to find YSOs associated with molecular 
clouds, so we start by removing stars. The foreground 
extinction, up to an allowed maximum Ay of 30 mag in 
this study is accounted for with an extinction law de- 
rived from GLIMPSE data (jlndebetouw et alj 120051 ) in 
the process of fitting stellar photosphere SEDs. We cal- 
culate a best-fit x 2 normalized by the number of data 
points between 4 and 8 used in the fit. We consider all 
sources with \ 2 /Ndata < 4 as good fits to reddened stel- 
lar photospheres. The total of sources with Ndata > 4 
separates into 69,316 well fit sources and 1,451 badly 
fit sources by this threshold value. In most cases, even 
highly reddened stars are classified as good fits because 
the fitting tool includes the extinction law. 
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TABLE 4 

Number of YSO candidates for Stages and Classes 



Region 



No. 



Physical 



Observational 







Stage 0/1 


11 


111 


Amb. 


Class 1 


Flat 


11 


ill 




Sources identified as 


YSOs with the SED fitter 








All^ 
/ill 


737 


228 


255 


5 


249 


256 


243 


Zoo 


r> 

o 


Fore. b 


300 


100 


93 


5 


102 


57 


75 


165 


3 


W51 c 


437 


128 


162 





147 


199 


168 


70 





Outer d 


281 


64 


137 





80 


88 


132 


61 





W51A ° 


72 


31 


9 





32 


52 


16 


4 





W51B f 


84 


33 


16 





35 


59 


20 


5 







Sources 


detected in 


all four 1HAC bands and MIPS 24 fim 






All 


305 


61 


196 


4 


44 


85 


108 


112 





Fore. 


109 


26 


61 


1 


18 


17 


23 


69 





W51 


196 


35 


135 





26 


68 


85 


43 





Outer 


164 


24 


118 





22 


45 


78 


41 





W51A 


15 


8 


6 





1 


12 


2 


1 





W51B 


17 


3 


11 





3 


11 


5 


1 







Sources classified as 


YSOc by Harvev et al. 


(2007) criteria 






All 


270 


59 


169 


1 


41 


83 


97 


90 





Fore. 


101 


26 


57 


1 


17 


17 


23 


61 





W51 


1(59 


33 


112 





21 


66 


71 


29 





Outer 


137 


22 


95 





20 


13 


67 


27 





W51A 


15 


8 


6 





1 


12 


2 


1 





W51B 


17 


3 


11 





3 


11 


5 


1 






a Entire sample of YSOs identified by SED fitter 
b Foreground YSOs with A v < 10 mag 
c W51 YSOs with A v > 10 mag 

YSOs with Av > 10 mag and located in the outside region excluding the central region 
(48.80° < I < 49.65° and -0.48° < b < -0.15°) 

c YSOs with A v > 10 mag around W51A (49.30° < I < 49.65° and -0.48° < 6 < 
-0.15°) 

f YSOs with A v > 10 mag around W5 IB (48.80° < I < 49.30° and -0.48° < b < -0.15°) 



Next, we fit 1,451 sources that are poorly fit by red- 
dened stellar photosp heres using the YSO models of 
iRobitaille et alJ (|2006[ L We allow the distance range to 
be from 5 to 9 kpc and the interstellar extinction from 
to 60 mag in V band for fitting YSOs. In the same man- 
ner as fitting stellar photospheres, we consider sources 
with x 2 / N data < 4 to be well-fit. This YSO fitting pro- 
cess separates these 1451 sources into 1327 good fit YSO 
sources and 124 bad fit sources. Among good fit sources 
we discard those not detected at 24 /jm but which show 8 
(im flux excesses above an extrapolation of the 3 shorter 
wavelength IRAC bands. These sources with 8 /im excess 
emission are usually found to be a noise peak or a diffuse 
background feature. In cases which show an IR excess in 
the MIPS 24 /im band only, we also discard the source, 
because these may be a false match due to the difference 
in resolution between GLIMPSE and MIPS GAL. We ex- 
tract 14 bright sources which may be evolved stars on the 
asymptotic giant branch (AGB). For example, carbon 
stars enshrouded by dusty envelopes have SEDs similar 
to YSOs. Most of the 124 objects classified as bad fits 
(% 2 1 N data > 4) are stars with saturated fluxes or signs of 
variability or bad photometry. After replacing the uncer- 
tainties of sources with questionable fluxes to upper or 
lower limits, we moved five YSO candidates among bad 
fit sources to the "well-fit" sample. Two samples were 
moved to the AGB candidate category. After all these 
adjustments, we distinguish 737 good YSO candidates, 
576 bad samples, and 16 bright sources which are likely 
to be AGB stars. 

We divide the final 737 YSO candidates into three 
evolutionary stages defined by the well-fit model disk 
mass Mdisk and the envelope accretion rate M env , both 



normalized by the mass of the central star M* ; Stage 
0/1 with Menv/M* > 1CT 6 yr" 1 , Stage II YSOs with 
Mcnv/M* < 1CT 6 yr" 1 and M disk /M* > 1CT 6 , and Stage 
III with M env /M* < 1(T 6 vr" 1 and M disk /M* < 1CT 6 . 
Stage 0/1, II, and III YSOs have significant infalling en- 
velopes, opti cally thick disks, and optically thin disks, 
respectively (jRobitaille et al.ll2006f h We determine the 
evolutionary stage of each source using the relative prob- 
ability distribution for the stages of all the "well-fit" 
models. The well-fit models of each source are defined 

by 



A 



A 2 - 

A mm 



< 2A^ dat a, 



(2) 



where Xmm ^ s ^ ne goodness-of-fit parameter for the best- 
fit model. The relative probability of each well-fit model 
is estimated according to 



P{X 2 ) = e-(* 2 - x » 



,)/2 



(3) 



and is normalized. After a probability distribution for 
the evolutionary stage of each source is constructed from 
the Stages of all the well- fit models, the most probable 
stage of each source is determined by the condition that 
SP(Stage) > 0.67. If this condition is not satisfied, then 
the stage of the source is classified as ambiguous, though 
the object is still counted as an YSO. 

From well-fit models for each source derived from the 
SED fitting tool, we calculated the x 2-we ighted model 
parameters: interstellar extinction to the source (Ay), 
stellar mass (M*), and total luminosity (L to t ■ total 
bolometric luminosity, including luminosity produced by 
accretion). Note that Ay is the foreground interstellar 
extinction, which does not include the extinction pro- 
duced by a circumstellar disk or envelope. The \ 2 ~ 
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Fig. 3. — Top: Example plots of rejected sources from final sam- 
ples. Left: Source with IR excess emission only in the IRAC 8.0 fim 
band, but not detected in the MIPS 24 (im band. Right: Source 
with contamination by diffuse background emission in 5.8 /^m and 
8.0 fim bands. Bottom: Example SEDs of AGB candidates. Tri- 
angles are datapoints with saturated fluxes that were used as lower 
limits for the model fitting. 

weighted average Ay distribution for all YSO candi- 
dates shows the linear variation in the direction of in- 
creasing extinction. We consider 437 YSOs with more 
than 10 mag Ay of foreground extinction as sources as- 
sociated with the W51 complex because the foreground 
extinction toward the W51 compl ex is Ay ~ 20 mag 
(|Goldader fc Wvnn-Williamj[l99l . The majority (97% 
(192/198)) of the YSO candidates detected in all 2MASS 
bands have Ay values < 10 mag, and are classified as 
foreground objects; while 95% (295/312) of the YSO can- 
didates detected only in IRAC/MIPS bands have fitted 
Ay > 10 mag and are classified as associated with the 
W51 region. We list magnitudes and classifications for 
the entire list of 737 YSO candidates in Table [2] . Fore- 
ground sources and W51 sources are denoted as "F" and 
"W" , respectively. We present some of the derived phys- 
ical parameters from the SED fitting for all 737 YSOs 
in Table [3] We report lower and upper ranges of 68% 
probability for each candidate YSO. In cases where the 
range in the 68% probability distribution is very small, 
we list ranges for the 95% confidence intervals and mark 
these with an asterisk. 

In the last three columns of Table 02 we give evolu- 
tionary stages and SED classes for each YSO candidates. 
We classify the 737 YSO candidates as 228 Stage 0/1, 255 
Stage II, 5 Stage III, and 249 ambiguous YSO candidates 
based on the disk mass Mdisk and the envelope accretion 
rate M env of each source. The spectral index chxrac is 
the value derived from the flux densities of the sources 
detected in all four IRAC bands. Among 737 YSO can- 
didates 561 sources are detected in all four IRAC bands 
and divided into 83 Class I, 164 Flat, 195 Class II, and 
19 Class III. Using the flux densities available between 2 
and 24 /im to derive the spectral index, ah sources 

are classified as 256 Class I, 243 Flat, and 235 Class II, 
and 3 Class III. We summarize the physical Stages and 
observational Classes of YSOs in Table [4] 

Figure EJJcf) shows the distribution of 561 YSO can- 
didates detected in all 4 IRAC bands in the IRAC 
[3.6] — [4.5] vs. [5.8] — [8.0] color-color diagram for com- 
parison with Figure [H 6) and (c). A total of 176 YSO 
candidates (24% of the total YSOs) are identified from a 



combination of JHK S and MIPS 24 /^m fluxes although 
these are not detected in all 4 IRAC bands. In Figure 
[2]Jcf), YSO candidates are classified according to their 
most probable evolutionary Stage: red for Stage 0/1, yel- 
low for Stage II, and cyan for Stage III. Since Stage III 
YSOs with optically thin remnant disks and their SEDs 
are dominated by photospheric emission, those objects 
are difficult to identify by IR excess alone. Despite a low 
probability, 5 Stage III sources are in our final samples 
because the evolutionary Stage is determined from mul- 
tiple YSO models of each YSO candidate statistically. 
Green dots are for ambiguous sources that are YSO can- 
didates but do not have a sufficiently well-defined prob- 
ability distribution to determine the evolutionary Stage. 
Blue dots are sources that were well-fit by YSO SEDs but 
were not included in the final sample. Figure [3l top shows 
examples of these rejected sources with IR excess emis- 
sion only in the IRAC 8.0 /im band (left) and with con- 
tamination by diffuse background emission (right). Ex- 
amples of SEDs for two AGB candidates are plotted in 
Figure [3l bottom. 

2.4. Estimating Contamination 

Because other IR-excess populations (e.g., AGB stars 
and galaxies) have colors or SEDs very similar to YSOs, 
the SED-fitting method will inevitably have contami- 
nation from these sources. Therefore we use an em- 
pirical method to estimate the contamination of YSO 
candidates from galaxies and dusty evolved AGB stars. 
GLIMPSE is a shallow survey of the Galactic plane while 
c2d is a deep survey of nearby star forming regions 
away from the Galactic plane: Ch ameleon II, Lupus , 
Perseus, Serpens, and Ophiuchus (|Evans et alj [2009*1 . 
The contamination by external galaxies is s erious in deep 
IRAC observations , such as the c2d Survey (jHarvev et alJ 
120071 : IPorras et al]l200l . In GLIMPSE, evolved stars 
like AGB stars with dusty winds and unresolved plane- 
tary nebulae (PNe) may be m ore serious contami nation 
sources than external galaxies (|Povich et al.ll2009h . 

To estimate an upper limit for evolved star (AGB) 
contamination, we use the color crit erion ([8.0] — 
[24] < 2.2 mag) of iWhitnev et alJ £2008). Their color- 
magnitude criterion derived from the LMC is not ap- 
propriate for W51, given the difference in distance and 
the fact that the AGB stars will not all be at the same 
distance as the YSOs. Note that relatively few extreme 
AGB stars are expected to populate the upper-right re- 
gion of the [8.0] vs. [8.0] — [24] color-magnitude diagram, 
and in the Galaxy such sources are likely to be saturated 
in GLIMPSE. Sources with [8.0] - [24] < 2.2 mag are 10 
of our 345 YSO candidates detected in [8.0] and [24] /im 
bands, or ~ 3% (Figure [4^6)). They do yield a similar 
(low) surf ace density of AGB stars as found for the M17 
region bv lPovich etaLl (|2009l ). The distribution of re d 
sources in the Galactic midplane (jRobitaille et al.l l 2008) , 
which consist mostly of YSOs and AGB stars, also shows 
that the number of YSO candidates represents a signifi- 
cant enhancement toward the W51 region while the sur- 
face density of candidate AGBs is very low in our region 
of 1.25° x 1.00° centered on (l,b) = (49.375°, -0.2°). 

For each source detected in all IRAC bands and 
in MIPS 24 /im, we calculated the probability of be- 
ing an extragalactic contaminant using the formulas in 
I Harvey et al.l (|2007l ). These formulas were derived for 
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Fig. 4. — (a), (ft), (c) Color-magnitude diagrams of YSO candi- 
dates in W51. The solid lines show the color-magni tude cuts for 
defini ng the YSO candidate criterion developed by IHarvev et al.l 
(2007) but with a total extinction of Ay = 25 mag along the line 
of sight, rather than the average value of 7.5 mag for the Ser- 
pens cloud. The dot-dashed lines show hard limits, fainter than 
which objects are excluded from the YSO category. Open sym- 
bols marked by circles and diamonds are objects which may be 
AGBs and extragalactic sources, (d) Plot of the number of sources 
vs. probability of being a background conta minant. The vertica l 
dashed line shows the separation chosen by Harvey et al. (2007) 
for YSOs vs. extragalactic candidates. 

the Serpens molecular cloud field, which has an average 
ex tinction of Ay ~ 7.5 mag, from the distribution found 
bv lHarvev et alj ([20071 ). The Serpens cloud lies at b = 5° 
latitude, so almost all of the extinction must be associ- 
ated with the cloud. In contrast, our W5I field is right 
in the plane of the Galaxy so the line of sight to any 
external galaxy should have a minimum extinction cor- 
responding to at least 25 kpc in the plane. We assume 
a minimum of 1 mag/kpc or Ay ^ 25 mag for our field 
or 17.5 mag more t han used for the Se rpens field. We 
therefore adjust the IHarvev et al.l (|2007l ) criteria by the 
expected additional extinc tion at [4.51, [8.0], and [24 1 , us- 
ing the extinct i on la ws of llndebetouw et al.l (|2005l ) and 
iFlahertv et all (pOOl . These give an additional 0.7 mag 
in all three of the IR bands. We add this amount to 
the constants for th e 3 IR magnitudes in the formulae of 
IHarvev et alj ([2007D . When we applied this modified set 
of color-magnitude criteria to 305 YSO candidates de- 
tected in all IRAC bands and in MIPS 24 /im, 35 sources 
were in the galaxy- candidate reg i on of the probability 
histogram using the IHarvev et al.l (|2007f ) cutoff of log P 
= —1.47. We note, however, that this cutoff represents 
a very conservative estimator, i.e., a probability as low 
as 0.034 counts as a galaxy candidate. Figure [4] presents 
three color-magnitude diagrams sh owing the pos s ibility 
of contamination using the modified IHarvev et alJ (|2007l ) 
criteria, where the solid lines from their formulae are 
shifted down by 0.7 mag in the magnitude axes. The 
joint probability histogram in Figure [UJd) shows that, of 
the 35 sources in the nominal GALc region, only 8 oc- 
cupy the bins with a probability greater than 0.06 (log 
P > —1.25) of being an extragalactic contaminant, i.e., 
about 3% of the sample of 305 objects. 

We therefore estimate an upper limit of about 6% to 
the contamination by galaxies and AGB stars in our list 
of YSOs. This limit implies a smaller contamination 



fraction compared to other regions, e.g., 10% for M17 
(|Povich et al.ll2009T ). 40% fo r the intrinsically red s ources 
in the Galactic midplane dRobitaille et ail [2008) , and 
55% for the LMC ([Whitney et al.ll2008ft . Therefore, con- 
sidering up to 2% contamination by unresolved Plane tary 
Nebulae (|Robitaille et alj|2008t [Whitney et al.ll2008D . we 
estimate an upper limit to the contamination of 8% in 
our YSO candidates. 

3. RESULTS AND DISCUSSION 

Our target area contains very active star-forming re- 
gions in the disk of the Milky Way (Figure [T]) . We find 
a total of 737 YSO candidates within the region around 
the W51 complex based on the SED fitter. After dividing 
737 YSOs into foreground and W51 YSOs based on the 
amount of foreground extinction, we focus on 437 YSO 
candidates associated with the W51 complex. 

These 437 sources divide into 199 Class I, 168 Flat, 
and 70 Class II when classified by the empirical SED 
spectral index, 0:2-24- For comparison, the percentage 
of sources in these observational SED classes for all c2d 
clouds is reported as 1 6% Class I, 12% F lat, 60% Class 
II, and 12% Class III ([Evans et al J 12009). The ratio of 
Class I plus Flat spectrum to Class II sources in W51 
is about 5, which is completely different from the low 
mass star-forming regions such as Cha II (0.16), Lupus 
(0.29), Perseus (0.57), S erpens (0.42), and Ophiuchus 
(0.47) (lEvans et all 120091) . The physical "Stages", I, II, 
and III should roughly corresp ond to the observat ional 
"Classes" I+Flat, II, and III ([Harvey et all 12001 . In 
W51 the 168 Flat sources divide into 30 Stage 0/1, 75 
Stage II, and 63 with ambiguous stage. Using the phys- 
ical evolutionary models for all sources associated with 
W51, 128 are classified as Stage 0/1, 162 as Stage II, and 
147 as ambiguous sources. The massive star-forming re- 
gion W51 appears to contain YSOs in a very early stage 
as indicated by both observational and physical classifi- 
cations. The ratio of Stage 0/1 to II for the entire set of 
YSOs in the large 1.25° x 1.00° region is similar to that 
for other massive star forming regions like M17 (25/40; 
i Povich et al]l2009h . N66 in the SMC (33/50:lSimon et all 
[2007D, and the LMC f 145 7147: IWhitnev et al.H2008fl . 

In the central active star-forming region (inside the 
boxed area in Figure the ratio of Stage 0/1 to II for all 
YSOs is 2.56 (64/25). In the same region this ratio for the 
massive YSOs (> 5M ) only is 2.71 (46/17). The ratio 
of YSOs outside the central bright region is 0.47 (64/137) 
and the ratio for sources greater than 5M Q in the same 
region is 0.92 (23/25). The increase in Stage 0/1 to II is 
due in part to a bias toward detecting sources with pos- 
itive spectral indices in regions of bright mid-IR nebular 
emission. M17 provides an extreme example of bright 
mid-IR nebular emission, where t he ratio reaches 1. 3 
(13/10) in the central bright region ([Povich et al.l [2009). 
In both W51 and M17, the relative fraction of sources 
with Ambiguous Stage determinations increases in the 
bright nebular regions, due to the decreased sensitivity 
of the MIPS 24 /im point-source detections, an impor- 
tant discriminant betwe en Stage 0/1 and Stage II sources 
([Indebetouw et al1l2007l ). However, these selection biases 
become less important for more luminous sources, and 
the large number of Stage 0/1 objects with M > 5 M Q 
found in our sample, 69, is unprecedented among studies 
of Galactic star formation regions to date. It implies that 
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Fig. 5. — (a) IRAC [3.6] — [4.5] vs. [4.5] - [8.0] color-color diagram with YSO candidates selected by the SED fitting tool. Points in 
red mark Stage 0/1, ye l low m ark Stage II, and green mark ambiguous sources. The solid line represents the region of the color selection 
criteria of ISimon et al.[ 120071). The dashed line sho ws the region shifted by the reddening vector for the general visual extinction toward 
W51 of 20 mag (Goldadcr & Wynn- Williams 1994). (b) IRAC [3.6] - [4.5] vs. [5.8] - [8.0] color-color diagram. The solid line represents 
the disk domain of Allen ct al. (2004). The dashed line shows the region shifted by t he reddening vector, (c) [3.6] — [5.8] vs. [8.0] — [24] 
color-color diagram. Reddening vector for A^=20 mag based on the extinction laws of Indcbctouw ct al. (2005) and Stutz et al. (2009, in 
preparation) is shown as an arrow. Black crosses show typical photometric errors. 



W51 is not only very massive but that it could also be 
exceptionally young, even when compared to M17. This 
result is consistent with the idea that W51 IRS 2E, the 
most luminous embedded source in W51A, is a highly 
embedded, very young O star that has yet to produce 
a det ectable hypercompact H II region (jFigueredo et al.1 
l200cl . 7 

3.1. Color- Color Diagrams 

Figure [5] shows the distribution of the final 437 YSO 
candidates in IRAC color-color diagrams. The sources 
plotted in Figure [5] are classified accordi ng to their most 
proba ble evolutionary Stage from the iRobitaille et alJ 
(2006) models. Red dots are Stage 0/1 YSOs that are 
relatively unevolved, heavily embedded in their natal en- 
velopes. Yellow dots are Stage II YSOs, like Class II 
T Tauri stars, with optically thick circumstellar disks. 
There are no Stage III YSOs with optically thin rem- 
nant disks in this sample. Green dots are objects with 
ambiguous stage classification. In Figure \Ma ), we see 
that the color selection criteria of lSimon et al.l (|2007[ ) in- 
cludes many reddened stars as well as YSO candidates 
within the solid line. If the color selection criteria are 
shifted by ~ 0.3 mag in the [3.6] — [4.5] color (correspond- 
ing to Ay = 20 mag), relatively clean YSO candidates 
are selected inside the das hed line. In Figure \5j b) , if the 
"disk domain" defined by lAllen et alJ (|2004D also moves 
vertically in the [3.6] — [4.5] color, many Stage II YSO 
candidates with optically thick circumstellar disks fall in 
the disk domain region. The shift of the color selection 
criteria in the direction of the reddening vector is con- 
sistent with the fact that the W51 complex is a more 
massive star forming region with m ore extinction tha n 
the star forming region studied by I Allen et al.l (2004). 
Figure [5]( c) shows the importance of the 24 /mi flux in 
classifying the evolutionary stages of YSO candidates. 
Stage 0/1 objects (red dots) are redder than Stage II ob- 
jects (yellow dots) in [8.0] — [24] color space. Most of 
the ambiguous sources (green dots) are in the boundary 
region between Stage 0/1 and Stage II. 

7 IRS 2E is saturated in the GLIMPSE and MIPSGAL images, 
hence it is not included in our sample of candidate YSOs. 



3.2. Spatial Distribution of YSOs 

Figure [6] shows the spatial distribution, classified by 
evolutionary stage, of the 437 YSO candidates associ- 
ated with the W51 region on the IRAC 8 fim image. 
YSO candidates are marked in red for Stage 0/1, yel- 
low for Stage II, and green for ambiguous sources by 
evolutionary stage. Stage 0/1 YSOs are clearly concen- 
trated toward the bright H II regions. Stage II YSOs 
are distributed outside of the central dense region. Most 
sources within the bright 8 /tm emission near the central 
region show strong clustering. These clustered sources 
are associated with the W51 molecular cloud. Several 
clusters of YSOs are also apparent in the outer region. 
A very strong concentration of YSOs is found in an IR 
dark cloud (IRDC) centered on (1, 6)=(49.4°, 0.0°). The 
detailed analysis of the CO observations and IRDC as- 
sociated with this YSO clustering will be presented else- 
where (Kang et al., in preparation). 

Figure [7] shows the spatial distribution of the YSOs 
classified by the mass of the central star on the IRAC 8 
/im image. The vicinity of the compact radio continuum 
sources marked as crosses in Figure [T] is bright in the 
IRAC 8 /im band due to PAH emission features stimu- 
lated by UV radiation from H II regions. The mass range 
of YSOs is from 1 to 18 M©. We divide all YSOs into 
three mass ranges: M* < 5M Q , 5M Q < M* < 8M Q , 
M* > 8M . Most of the massive YSOs (> 8M Q ) are 
located near the main H II regions of W51. Cluster- 
ing of massive YSOs is clearly apparent in the central 
region, while the YSO cluster centered on (1, 6)=(49.4°, 
0.0°) consists of many lower mass YSOs and only one 
massive YSO. The mass distribution of YSOs shows evi- 
dence of mass segregation with the most massive objects 
exhibiting a strong concentration close to the center of 
the H II regions. We can quantify this trend by comput- 
ing the mass function of YSOs for various regions within 
our target area. 

3.3. Mass Function of YSOs 

We investigate the variation with spatial distribution of 
the mass function (MF) of YSOs. We divide our target 
field (1.0° x 1.25°) into two regions: the central part 
including very bright emission in all bands from near-IR 
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Fig. 6. — Spatial distribution of YSOs in W51, overplotted on IRAC 8.0 £im image. YSO candidates are marked in red for Stage 0/1, 
yellow for Stage II, and green for ambiguous sources. 
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Fig. 7. — The mass distribution of 437 YSOs associated with the W51 region overplotted on the IRAC 8.0 £tm image (This is the same 
image as in Figure \E\ but displayed with inverted grayscale image to show more clearly the structure of very strong PAH emission in black). 
Small (cyan), medium (yellow) and large (red) dots represent YSOs with M < 5A/q, 5Mq < M < 8Mq and M > 8Mq. Box shows the 
boundary of the central region for comparing the shape of the IMF. Left side of the box is the W51A region and right side is the W51B 
region. 
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Fig. 8. — The mass-luminosity relation for YSOs in the W51 
region. Blue and red dots represent foreground and W51 sources, 
respectively. Open circles are the sources used for MF de rivation. 
Green line shows the ZAMS mass-luminosity relation by Malkov 
H2007D for intermediate mass stars (1 — 30Mq). 

to 21 cm continuum (48.80° < I < 49.65° and -0.48° < 
b < —0.15° shown by the solid box in Figure [7]) and the 
outside region excluding the central region. We also find 
the MF of the W51A and W51B regions of the central 
part, separately. 

To obtain as clean a sample as possible, we remove 
from the final source list those objects with a probability 
of being extragalactic and AGB contaminants, as dis- 
cussed in Section l2~4l To compare MFs without bias of 
source selection we set two criteria. The sensitivity limit 
of IRAC sources is affected by the brightness of the dif- 
fuse emission, in particular at IRAC 8.0 /im. Therefore, 
the sensitivity limit in star-forming regions with bright 
diffuse emission is certainly d ifferent from the outside 
region. iRobitaille et all (|2008ft found that sources with 
lower limits on the IRAC 4.5 and 8.0 /im fluxes of 0.5 
and 10 mjy respectively are not affected by the varia- 
tions in the diffuse emission on a Galaxy-wide scale. We 
apply these lower limits on the IRAC 4.5 and 8.0 /xm 
fluxes to our final YSO sources, which selects 113 YSOs. 
The mass distribution of 113 YSOs turns over at 5 M . 
Therefore we use only those 94 YSOs with M* > 5 M 
for comparing the shape of the mass function of YSOs 
in the W51 region to avoid introducing bias due to sen- 
sitivity and completeness. Although we have used very 
conservative brightness cutoffs in constructing our mass 
functions, the mass functions are unlikely to be 100% 
complete in the central regions, due to the selection ef- 
fects imposed by confusion and extremely bright nebular 
emission. In the case of more luminous objects, the stel- 
lar luminosity dominates in the total luminosity of the 
source, while the accretion luminosity is clearly negligi- 
ble, so the SED-determined luminosity simply sets the 
corresponding stellar mass given the distance. We find 
that the mass-luminosity relation for the sources used in 
the MF derivation agrees with a recent re-evalu a tion o f 
the ZAMS mass-luminosity relation by iMalkovl (|2007l) . 
as shown in Figure[8l The mean distance of 94 YSOs de- 
rived from the SED fitter is 6.6 ± 0.2 kpc which is close 



to the 6.1 ± 1. 3 kpc measured for water masers in the 
W51A region (jlmai et al.H2002T ). 

The mass distribution from 94 YSOs which satisfy the 
sensitivity and completeness criteria in our target region 
has a best-fit slope of the MF, T = -1.70 ± 0.10, as 
defined by 

d£(M*) 



Ml oc 



d(logAf*) 



for Mi, > 5M Q (solid line in Figure [5K a) ) . This value 
is very c lose to the st ellar IMF slope of T = —1.8, de- 
rived by IRanal (|1987f ) for the solar neighborhood field 
stars (1.5Mq — 100M Q ). Foreground sources which sat- 
isfy the same sensitivity and completeness have the MF, 
T = -2.13 ± 0.29 (dash-dotted line in Figure Efa)). 
However the YSO MFs of the central active star form- 
ing region (r = —1.26 ±0.12) and the outer region 
(r = —2.36 ± 0.26) are significantly different in Figure 
\§lb). The MF for the inner region has a slope similar to 
the Sal peter law, while the MF of the outer part is very 
steep. iKroupa fc~W cidncr (|2003l ) note that field-star 
IMFs are always steeper for masses over 1M than the 
stellar IMF derived from an individual massive s tar clus- 
ter. In the Large and Small Magellanic Clouds, IMassevi 
(|2002f) found a very steep IMF slope (r 4 ± 0.5) out- 
side of the OB associations. The IMF is often observed 
to be shallower in the central regions of clusters. This 
mass segregation effect might occur because the most 
massive stars are either born near the center or migrate 
toward the center. The question is then whether this 
segregation is due to initial conditions or to dynamical 
evolutionary effects. In the W51 H II region complex, 
most of the massive YSOs in the early stages of evolu- 
tion are concentrated in the central region (see Figure [6] 
and [7J. Because we derive the MFs using YSOs in the 
early stages, still surrounded by disks and/or infalling en- 
velopes, these YSOs must be located in their birthplaces. 
Therefore the spatial variation of the MF in W51 indi- 
cates that the mass segregation is caused by initial condi- 
tions, not dynamical effects. The difference between the 
central and outer regions is whether the sources are as- 
sociated with the most active star-forming region or not. 
In the outer region YSOs may be forming in relatively 
small or isolated molecular clouds. 

We compare MF of YSOs calculated from Spitzer ob- 
servation to the stellar mass function. IMFs of many 
dense clusters in the solar neighborhood are close to 
the Salpeter IM F slope, T = — 1.3 5 on a plot with log- 
mass intervals. iScalol (|1986l |1998[ ) reviewed the deter- 
mination of the stellar IMF extensively and concluded 
the IMF slope of field stars to be -1.7 ± 0.5 in the 
range of 2 - 1OM . INinkov et all (|1995D found a slope 
of T = -1.38 ± 0.19 for masses between 2.5 - 3OM 
from the OB cluster in the H II region IC 1805 in the 
Cas OB6 association associated with the molecular cloud 
W4. They noted that the mass function determined from 
young galactic OB clus ters varies from T = —1.0 to —1.4. 
IGarmanv et al.l (|1982f) determined a slope of the IMF 
of r = —1.6 for massive stars (> 2OM ) within 2 kpc 
of the Sun. They found evidence that the slope of the 
IMF varied with galactic radius: the IMF of stars inside 
the solar circle is T = —1.3 and outside is T — —2.1. 
They suggested that the excess of massive stars inside 
the solar circle was due to the OB associations in the 
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Fig. 9. — (a) The mass function for YSO candidates after selecting only > 0.5 and 10 mjy at 4.5 and 8.0 and cut off > 5 Mq. The 
solid line and the dash-dotted line represent the mass function of YSOs associated with W51 and foreground YSOs, respectively. (6) The 
mass function of the central region (box in Figure [7} (solid line) and the outside region (dash-dotted line), (c) The mass function of the 
W51A (solid line) and the W51B (dash-dotted line) region. All histograms are shown with y/N errorbars. 



Carina and Cygnus spiral arms. iHunter et alj (|1996l ) es- 
timated an IMF slope of — 1.6±0.7 for intermediate-mass 
stars (6.5 - 18 My ) from NGC 604 in the galaxy M33. 
iMassev fc Hunter! (|1998l ) derived a slope close to the 
Salpeter value from the 30 Doradus region of the LMC, 
suggesting that ther e is no strong depen dence of the IMF 
slope on metallicity. iSabbi et al.l (|2008f ) found a present- 
day mass function of T = -1.43 ± 0.18 of NGC 346, the 
active star-forming region in the SMC. lHomeier fc Alvesl 
(|2005f ) derived a slope of r = — 1.6 ± . 3 in t he giant ra- 
dio H II region W49A. lOkumura et al.l (|2000l) found that 
the IMF of the massive stars (M > 4M Q ) in G49.5-0.4 of 
W51A is consistent with an IMF slope of —1.8, but that 
the IMF shows a statistically significant excess of stars 
in the high mass range above 30M©. These various ob- 
servations suggest that low density regions have slightly 
steeper IMFs than the Salpeter slope, while more mas- 
sive clusters have a mo re massive upper end of the IMF 
( Elmegreen 2004, 2006) . The results we present here ap- 
pear to support this conclusion. 

We divide the central active star-forming region into 
W51A and W51B (Figure and find the slopes of the 
MF of T = -1.17 ± 0.26 and T = -1.32 ± 0.26, re- 
spectively (Figure Hfc)). The slopes are the same as 
the Salpeter IMF within the errors. Similar results that 
the MFs of YSOs are consistent with the Salpeter IMF 
are reported from o ther active star-for ming regions, e.g., 
M17 in the Galaxy (|Povich et a l.l l2009ri. N66 in the SMC 
(|Simon et al.ll2007h . and the LMC (|Whitnev et al.ll2008l) . 
It is not surprising that our MF of YSOs in this active 
star-forming region is similar to the stellar IMF because 
even the core mass function, which has T = —1.3 for 
M > 0.8M© in Perseus, Serpen s, and Ophiuchus, is con- 
sistent with the Salpeter IMF (jEnoch et al.ll2008f ). In a 
future paper, we will present the detailed properties of 
molecular clouds associated with the YSOs in W51A and 
W51B using our 12 CO and 13 CO J = 2 — 1 emission line 
maps of the region studied in the present work. 

4. SUMMARY 

We have found 737 candidate YSOs in a 1.25° x 1.00° 
area near the W51 GMC using Spitzer Space Telescope 
data. We distinguish YSOs from main-sequence or red gi- 



ant stars, asymptotic giant branch stars, and background 
galaxies by fitting model SEDs to fluxes of sources. We 
divide the total of 737 YSOs into two groups based on 
interstellar extinction for each YSO: sources associated 
with the W51 region and foreground sources. We iden- 
tify 69 highly embedded Stage 0/1 candidate YSOs in 
our field with masses > 5 M© (corresponding to mid- to 
early-B main-sequence spectral types) 46 of which are lo- 
cated in the central active star forming regions of W51A 
and W51B. From the YSOs associated with W51, we 
find evidence for mass segregation showing that the most 
massive YSOs are concentrated on the W51 H II region 
complex. We find a variation in the spatial distribution of 
the MF of YSOs in the mass range between 5 Mq and 18 
Mq . The mass distribution from YSOs in our entire tar- 
get region has a slope of T = — 1.70±0.10. Whenthesam- 
ple is divided into the active star-forming region and the 
surrounding outer region, the slope of MFs in the star- 
forming region (r = —1.26 ±0.12) is shallower than that 
of outer region (r = —2.36 ± 0.26). However, we do not 
find any difference in slopes between two active star form- 
ing regions within the complex: T = —1.17 ± 0.26 and 
T = -1.32 ±0.26 in W51A and W51B, respectively. The 
variation of the MF for YSOs embedded in the molecu- 
lar cloud implies that the distribution of stellar masses 
in clusters depends on the local conditions in the parent 
molecular cloud. 
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